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Abstract: After a period of inflationary expansion, the Universe reheated and reached full ther- 
mal equilibrium at the reheating temperature Tr. In this work we point out that, in the context of 
effective low-energy supersymmetric models, LHC measurements may allow one to determine Tr 
as a function of the mass of the dark matter particle assumed to be either an axino or a gravitino. 
An upper bound on their mass and on Tr may also be derived. 
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1. Introduction 

Dark matter (DM) remains an unknown component of the Universe. While it has so far es- 
caped detection, its existence has been convincingly inferred from gravitational effects that 
it imparts on visible matter through rotational curves of spiral galaxies, gravitational lens- 
ing, etc, [|l|. The effects of dark matter also can be seen on large structure formation and on 
anisotropy of the cosmic microwave background (CMB). The CMB in particular provides 
a powerful tool for determining the global abundance of cold DM (CDM). Recently, the 
Wilkinson Microwave Anisotropy Probe (WMAP) has performed a high-accuracy measure- 
ment of several cosmological parameters. In particular, the relic density of CDM has been 
determined to lie in the range Q 

17cdm/i^ = 0.104 ±0.009. (1.1) 

Since DM has to be electrically and (preferably) color-charge neutral, from the particle 
physics point of view, a natural candidate for DM is some weakly interacting massive 
particle (WIMP). Within standard cosmology, the WIMP is produced via a usual freeze- 
out mechanism from an expanding plasma. 

Among specific, well-motivated particle candidates for the WIMP, the by far most pop- 
ular choice is a stable lightest neutralino x of effective low-energy supersymmetry (SUSY) 
models. Most efforts have gone to exploring the lightest neutralino x of the Minimal Super- 
symmetric Standard Model (MSSM) as the lightest supersymmetric particle (LSP) that, 
in the presence of R-parity, is stable and makes up all, or most of, the CDM in the Uni- 
verse. In addition to an impressive experimental effort of direct and indirect searches for 
cosmic WIMPs, the Large Hadron Collider (LHC) will soon start exploring the TeV en- 
ergy scale and is expected to find several superpartners and to determine their properties. 
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In particular, some authors have explored the feasibility of determining the neutralino's 
relic abundance ^x^'^ from LHC measurements [|3|, Q]. Their conclusion was that, under 
favorable circumstances, this should be possible with rather good accuracy, of order 10% 
or better, although this may be challenging An analogous study has also been done 
in the context of the Linear Collider, where accuracy of a similar determination would be 
much better 

If a WIMP signal is detected in one or more DM detection experiments, and also at 
the LHC a large missing-mass and missing-energy signature, characteristic of the stable 
neutralino, is measured and implies a similar mass, and if perhaps additionally eventually 
its relic abundance is determined from LHC data, even if with limited precision, and agrees 
with the "WMAP range" (|l . l| ) , then the DM problem will most likely be declared solved, 
and for a good reason. 

However, such an optimistic outcome is by no means guaranteed. One realistic possibil- 
ity is that the neutralino, even if it is found as an apparently stable state in LHC detectors, 
may not be the true LSP and therefore DM in the Universe. Instead, it could decay in the 
early Universe into an even lighter, and (possibly much) more weakly interacting, state, 
the real LSP, outside the MSSM (or some other low-energy SUSY model) spectrum. In this 
case, current cosmic WIMP searches will prove futile, even after improving the upper limit 
on the spin-independent interaction cross section on a free proton 0"^^ from the current 
sensitivity of ~ 10~^ pb down to ~ 10"^'' pb, which is as far down as experiment can 
probably go given background from natural radioactivity. 

Moreover, the neutralino relic abundance, as determined at the LHC, may come out 
convincingly outside the range ( |1.1| ). In fact, a value of Q^^^ above about 0.1 can be easily 
explained in terms of a lighter LSP into which the neutralino, after its freeze-out, decayed 
in the early Universe Q. This is because the relic abundance of the true LSP is in this case 
related to fi^/i^ by the ratio of the LSP mass rriLSP to nix, "^hich is less than one. If ^x^'^ 



comes out below {'LI), several solutions have been suggested which invoke non-standard 
cosmology, e.g. quintessence-driven kination while preserving the neutralino as the 
DM in the Universe. However, if at the same time DM searches bring null results, this will 
provide a strong indication against the neutralino nature of DM. In contrast to the above 
attempts, we will consider a whole range of possible values of ^x^'^ LHC, both below 

and above 0.1. We will work instead within the framework of standard cosmology but will 
not assume the neutralino to be the DM in the Universe. 

In this context we point an intriguing possibility of determining at the LHC the re- 
heating temperature of the Universe. The framework we consider will therefore give us 
an opportunity to probe some crucial features of the early period of the Universe's his- 
tory. The reheating temperature Tr is normally thought of as the temperature at which, 
after a period of rapid inflationary expansion, the Universe reheated (or, more properly, 
defrosted), and the expanding plasma reached full thermal equilibrium. Determining Tr 
cannot be done with the neutralino as DM since it freezes out at Tf ~ m^/24, which is 
normally thought to be much below Tr.^ Here instead we assume a different candidate 



""^The possibility of Tr < Tt has been explored in and more recently in In this case one could 
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for the LSP (assuming R-parity) and cold DM, whose rehc density depends, at least in 
part, on Tr. This is the case for either an axino or a gravitino. The spin-1/2 axino (the 
fermionic superpartner of an axion) and the spin-3/2 gravitino (the fermionic superpartner 
of a graviton) are both well-motivated. The former arises in SUSY extensions of models 
incorporating the Peccei-Quinn solution to the strong CP problem. The latter is an in- 
herent ingredient of the particle spectrum of supergravity models. Both, like the axion, 
form a subclass of extremely weakly interacting massive particles (E-WIMPs) |12|. The 



characteristic strength of their interactions with ordinary matter is strongly suppressed 
by a large mass scale, the Peccei-Quinn scale fa ~ 10^^ GeV in the case of axinos and 
the (reduced) Planck scale Mp ~ 2.4 x 10^^ GeV for gravitinos. The mass of the axino 
rria is strongly model-dependent and can take values ranging from keV up to TeV |l3|| . 
The mass of the gravitino in gravity-mediated SUSY breaking schemes is given by 
niQ ~ Mg/Mp, where Ms ~ 10^^ GeV is the scale of local SUSY-breaking in the hidden 
sector, and is expected in the GeV to TeV regime. In other schemes of SUSY breaking mg 
can be (much) smaller. In this work we want to remain as model-independent as possible 
and will treat and as free parameters. Both axinos and gravitinos can be produced 
in decays of the neutralino (or another ordinary superpartner, e.g. the stau) after freeze- 
out, as mentioned above, or in thermal scatterings and decay processes of ordinary particles 
and sparticles in hot plasma at high enough reheating temperatures Tr - hence their relic 



density dependence on Tp. The possibility of axinos in a KSVZ axion framework |14| as 



cold DM was pointed out in ||8|, |T^] and next studied in several papers |16, |T^, 18, 19 1 



while axinos as warm DM was considered in . The gravitino as a cosmological relic was 
extensively studied in the literature, starting from |21, 22|, more recent papers on gravitino 



CDM include [E3L 24, 25, M, 27 



For definiteness, in this work we will first assume the lightest neutralino to be the 
lightest ordinary superpartner and the next-to-lightest particle (NLSP), although below 
we will also consider the case of the stau. Our main result is that, assuming the axino 
or gravitino as the true LSP and CDM, and that the (apparently stable) neutralino is 
discovered at the LHC and its "relic abundance" Q^h"^ is determined from LHC data, then 
one should be able to determine the reheating temperature in the Universe as a function 
of the LSP mass, or at least place an upper bound on it, as we show below. In the regime 
where thermal production dominates, in the axino case we find Tp cx fa/f^a while in the 
gravitino case Tp oc mg. Alternatively, in the non thermal production dominated regime, 
we find an upper bound on the allowed mass range of a CDM particle. Furthermore, fi-^/i^ 
at the LHC can be expected to come out either below or above (or for that matter even 
accidentally agree with) the "WMAP range" (|1.1|). 

The same holds true also in the even more striking case of the lighter stau taking the 
role of the NLSP instead. In this case the very discovery of an (apparently stable) charged 
massive particle at the LHC will immediately imply that DM is made up of some state 
outside the usual spectrum of low-energy superpartners. In this case, dedicated studies of 
the differential photon spectrum may allow one to distinguish between the axino and the 



think of determining experimentally Tr even with the neutralino as DM. 
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gravitino LSP 

We stress that, a detection of a (seemingly) stable neutralino at the LHC will not be 
sufficient to prove that the neutralino is the LSP and the CDM.^ Even establishing an 
apparent agreement of Vl^h? , as to be determined at the LHC, with the range will 
not necessarily imply the neutralino nature of DM, although admittedly it will be very 
persuasive. For this, a signal in DM searches must also be detected and the resulting 
WIMP mass must be consistent with LHC measurements of the neutralino mass. 

The paper is organized as follows. In section ^, we review the thermal and non-thermal 
production of axinos and gravitinos and next explain our strategy for determining at 
the LHC. In sections ^ and ^, we discuss the reheating temperature determination with 
axino LSP and gravitino LSP respectively, assuming the neutralino to be the NLSP. In 
section ^ we consider instead the lighter stau as the NLSP. We make final remarks and 
summarize our findings in section |6[ 



2. Cosmological production of relic axino or gravitino E-WIMPs 

First we briefly review main mechanisms of producing relic axino or gravitino E-WIMPs, 
assumed to be the LSP, in the early Universe. Since both are neutral Majorana particles, 
they are produced in an analogous way. Since the interactions of axino and gravitino 
are strongly suppressed with respect to the Standard Model interaction strengths by fa 
and Mp, respectively, they can be in thermal equilibrium only at very high temperatures 
of order 10^^ GeV. Axinos and/or gravitinos produced at lower temperatures are out-of- 
equilibrium, thus their abundance depends on the reheating temperature. 

Neglecting the initial equilibrium population (which we assume to have been diluted 
through inflation), we consider here two generic ways to re-populate the Universe with 
axinos or gravitinos. One proceeds via scatterings and decay processes of ordinary particles 
and sparticles in thermal bath. Its efficiency is proportional to their density in the plasma 
which is a function of Tr. Following we call it thermal production (TP). The other, 
dubbed non-thermal production (NTP) refers to (out-of-equilibrium) decays of the NLSPs, 
after their freeze-out, to E-WIMPs. In addition, there could be other possible mechanisms 
contributing to E-WIMP population, e.g. through inflaton decay, but they are much more 
model dependent and will not be considered here. 



Thermal production at high Tr. In considering thermal production of E-WIMPs 
we compute their yield ^^lsp = "-lsp/s, where ulsp is their number density and s is the 
entropy density. To this end, we integrate the Boltzmann equation up to the reheating 
temperature Tr. Their thermal abundance is then given by 

QJsV = mLSPlS"p^, (2.1) 



^Note that, for the neutrahno, or any other superpartner, to appear stable in an LHC detector, it is 
sufiicient that its Ufetime is longer than a microsecond or so. Of course, if it is unstable, then it will be 
immediately clear that the neutralino is not cosmologically relevant. 
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Figure 1: The thermal yield of axinos Y^^ (left panel) and gravitinos (right panel) as a 

function of Tr. In both panels the neutralino mass is — 300 GeV, the stau mass is m-f-^ — 
474 GeV, while gluino, squark masses are set at 1 TeV. In the left panel we take fa = 10^^ GeV; 
the yield scales as 1//^. The black solid curve has been obtained within the effective (plasmon) 
mass approximation, as implemented in fill , which we use in this work. For comparison, the blue 
curve has been obtained in by applying a hard thermal loop (HTL) technique. We also show 
the relative contributions to from squark, gluino and neutralino decays. In the right panel 
we compute the gravitino yield following ref. [ p5[ where the HTL technique was applied. In both 
panels the results obtained within the HTL formalism are more correct at Tr > 10^ GeV (dashed) 
but less applicable at lower Tr (dotted). 



where h is the dimensionless Hubble parameter and pcrit is the critical density. In other 
words, mLSP^i^p ~ '^•'^ ^ 10"^'^ (J7lsp^^/0-1) ' which is actually true in general. 

In the case of axinos we follow the procedure described in detail in |^ |l5|, |l^, |l^] . The 
main production channels are the scatterings of strongly interacting (s)particles described 
by the dimension-five axino-gluino-gluon term 

/:a^^Y^575[7^71rC• (2.2) 

The most important contributions come from 2-body processes into final states, « + j — > 
a + • • • [l^. Axinos also can be produced through the decays of heavier superpartners 
in thermal plasma. In this work we include gluino, squark, slepton and neutralino decays. 
At low reheating temperatures, comparable to the mass of the decaying particle, axino 
production from squark decay in thermal bath gives dominant contribution [p!5| ]. 

In the gravitino case, the analogous dominant goldstino-gluino-gluon dimension-five 
term is given by |23, S^] 



where the goldstino ipc is the spin-1/2 component of the gravitino. Note the dependence 
on mg/rriQ, while no analogous term appears in the axino case. In computing the gravitino 
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yield from thermal processes we follow ref. [25|. On the other hand, we do not include 



gravitino production from sparticle decays in the plasma |19| which change the yield by a 
factor of about two. 

Fig. IT] shows the yield of axinos (left panel) and of gravitinos (right panel) 
from thermal production as a function of Tr . Input parameters are given in the caption. 
As expected, in both cases the yield grows with increasing Tr. In the axino case, a sharp 
drop-off below Tr ~ 1 TeV is due to Boltzmann suppression factor exp(— m/T), with m 
denoting here squark and gluino mass; at lower Tr superpartner decay processes become 



dominant but are less efficient [15|. 



In both axino and gravitino production from scatterings in thermal plasma, some 
classes of diagrams suffer from infra-red divergence. An early remedy in terms of a plasmon 



mass as a infra-red regulator in the gravitino case in |22| was more recently improved using 
hard thermal loop (HTL) technique and a finite result was obtained |2j, The HTL 
technique is fully applicable only in cases corresponding to rather high Tr > 10^ GeV for 
which, in the right panel of fig. |^, the gravitino yield is marked with a dashed line for a few 
representative choices of mg. At lower Tr we use the same formalism but the result is less 
reliable (dotted parts). In the axino case an analogous calculation was performed in |18| 



but is of more limited use since for the axino to be CDM one requires Tr < 10 GeV |15] . In 



this work in computing the axino yield we therefore follow [15| and use the plasmon mass as 



a regulator. The result is shown as solid line in the left panel of fig. Q and labelled "effective 



mass approximation". For comparison, the blue line shows the result obtained in |18| using 
the HTL technique which is more correct than our treatment in the regime marked by a 
dashed line but less reliable in the one marked by a dotted one (Tr < 10^ GeV). In any 
case, the difference is of order a few which is not important for our purpose. 

Since the axino yield from thermal production is basically independent of the 
axino mass, until it becomes comparable to the masses of MSSM sparticles, its relic abun- 
dance ^J^h"^ is proportional to rria and can then be expressed as [^] 

■>TP,2 _ n 1 /' ^ f 



At high enough Tr thermal production is dominated by scattering processes and an appli- 
cation of the HTL leads to the following formula 



gs J vO.lGeV/ \ h ) V 10^ GeV 

where gs is temperature-dependent strong coupling constant, which in the above expression 
is evaluated at Tr. Note that, Y"^^ oc T-^/ f'^, as expected. 

In the gravitino case, due to the above-mentioned gravitino mass dependence in the 
denominator of the dimension-five terms of the gravitino Lagrangian, Y~^ oc l/m^'^, and 



the relic abundance calculated using a HTL technique can be expressed as |25] 



ni^h^ ~ 27 (-^] r ^QQ^^v ^ f (2 6) 
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where mg{fi) stands for the gluino mass evaluated at a scale fi ~ ITeV. 

Non-thermal production from NLSP decays. As the Universe cools down, all heav- 
ier superparticles first cascade decay into the lightest superpartners of the low-energy SUSY 
spectrum, in our case the NLSPs. The NLSPs then freeze out from thermal equilibrium 
and only later decay to LSP axinos or gravitinos. Since all the NLSPs decay into axino or 
gravitino LSPs, their number density is the same as that of decaying NLSP, Ilsp = ^nlsp, 
and the LSP relic abundance from non-thermal production is given by 

^LSp"^^ = ^^^f^NLSP^^- (2-7) 

"T-NLSP 

Once super symmetric particles are found at the LHC and relevant parameters of the 
NLSP are determined, then its freeze-out relic abundance can be calculated by solving the 
Boltzmann equation (under the assumption that it is actually stable until after its freeze- 
out).^ Since the NLSPs do not constitute CDM, we treat the quantity I^nlsp^'^ as a free 
parameter which will be determined experimentally at the LHC. 

Relic density and the LHC. The total abundance of the LSPs is the sum of both 
thermal and non-thermal production contributions 

n^sph' = nliph' + n^i^h\ (2.8) 

Since it is natural to expect that the LSP makes up most of CDM in the Universe, we can 
re-write the above as 

^LSph^ (7R,mLSP,mg,mNLSP, • • •) + ^^^^ ^NLSP^^ = ^^LSP^^ = ^CDuh^ ^ 0.1. (2.9) 

"INLSP 

This is our master formula. Once the neutralino NLSP is discovered and its mass is 
determined at the LHC with some precision, and so also J^nlsp^^ = f^x^^' then eq. ( p^ ) 
will provide a relation between Tr, and m-LSP- 

More specifically, in order to evaluate ^x^"^ the LHC, the parameters determining 
the neutralino mass matrix will have to be evaluated through measurements at the LHC, 
as well as the masses of some other superpartners and Higgs boson(s) Q. For our 
purpose, also the mass of the gluino rrig will have to be known since it determines the 
efficiency of E-WIMP production in thermal scatterings. We assume the usual gaugino 
mass unification which, at the electroweak scale implies Mi = 5/3tan^ 0vK-^2 = oti/asrrig 
(or Ml ~ O.5M2 — mg/6.5), where Mi and M2 are the bino and wino mass parameters, 
respectively. (If the gluino mass is actually measured, we can use it in our expressions 
instead.) In this work we assume the lightest neutralino to be mostly a bino in which case 

— Ml ~ mg/6.5. 

Note that, if I^nlsp^^ is determined to be larger than 0.1, the relic abundance of the 
LSP may easily agree with ( |l.lD by taking an appropriate mass ratio m-LSp/^iNLSP and 



Obviously, in our case f^NLSpft is not the true cosmological relic abundance of the NLSPs. Instead it 
is simply given by JInlsp/i^ = TiNLSpyNLSP^p^ where Fnlsp is the NLSP yield at its freeze-out. Here, for 
convenience we will keep calling JInlsp^^ a "relic abundance". 
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^NLSph^ m, (GeV) 

Figure 2: Left panel: Tr vs. ^nlsp^^ for mNLSP — 300 GeV and for ma — 0.01 GeV (solid blue) 
and ma = 1 GeV (dashed red) . The bands correspond to the upper and lower limits of dark matter 
density in (^Tl|). Right panel: Tr vs. ma for f^NLSP^^ — 100 (dashed blue), 0.1 (solid red) and 0.01 
(dotted black) . To the right of the solid vertical line the axino is no longer the LSP. In both panels 
we set fa 10" GeV. 

small enough Tr in order to suppress TP. Should, on the other hand, ^x^"^ came out to 
be smaller than 0.1, a substantial contribution of LSP production in thermal processes at 
large enough Tr will be sufficient. Actually, TP may be dominant even if Q^h'^ is larger 
than, or for that matter is close to, 0.1, depending on the LSP and the neutralino NLSP 
mass ratio. 

3. Axino dark matter 

We first consider the axino as the LSP and the dominant component of DM in the Universe. 
As mentioned earlier, for definiteness we take the NLSP to be the lightest neutralino; the 
stau case will be considered below. It is reasonable to expect that LHC experiments will 
be able to probe neutralino (and stau) mass ranges up to some 400 GeV, depending on 
other SUSY parameters. Here we are not concerned with experimental uncertainties of 
LHC measurements but rather illustrate the principle of estimating Tr for a given DM 
candidate in terms of relation ( p.9[ ). 

In the left panel of fig. ^ we show ranges of f^NLSP^^ ^-nd Tr, such that ila^^ is in the 
range ( |1 . l\ j . We take a fixed neutralino NLSP mass jtinlsp = 300 GeV and rria = 0.01 GeV 
(solid blue) and 1 GeV (dashed red). For each choice of the axino mass the two lines 
correspond to the upper and lower limits of dark matter density in (|1.1|). On the left 
side non-thermal production is subdominant while thermal production of axinos gives the 
main contribution to their relic density. In this case Tr corresponds to the value for which 
JlJ^/i^ ~ r^cDM^^, and depends on the axino and NLSP masses since Y^^ oc T^/ f"^. As 
Q^'^^h? increases, Vt^^h? must decrease (in order for the total abundance to remain close to 
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Figure 3: Maximum values of mLSP as a function of J^nlsp^^ for representative values of rriNLSP- 
Once both r^NLSP^^ and jtinlsp are determined from experiment, the upper bound on jtilsp can 
be derived. The plot applies both to the axino and to the gravitino LSP. 



0.1) and at some point becomes subdominant. This point is marked by an abrupt turnover 
of the contours from horizontal to vertical. 

In the right panel of fig. |2| we plot rria vs. Tr for the same NLSP mass and for 
^^NLSP^^ = 100 (dashed blue), 0.1(solid red) and 0.01 (dotted black). On the left side, 
where rria is small, TP dominates, Oj^/i^ ~ r^cDM^^) hence we find 

oc fl/rua. (3.1) 

In this regime is inversely proportional to the axino mass (for a given fa, see below). 



(This dependence can be seen analytically by plugging eq. (|2.5| ) into eq. (2^) and taking the 



limit of negligible contribution from NTP.) The relation (3T) allows one to derive an upper 
hound on Tr if we use the fact that axinos have to be heavy enough in order to constitute 
CDM. Assuming conservatively that > 100 keV |5| we find T^''^ < 4.9 x 10^ GeV. 

At larger rria the NTP contribution becomes dominant and the dependence on Tr 
is lost but one can still place a lower bound on Tr. Since in this regime the LSP mass 
becomes largest, this allows one to derive an upper bound on nia- This is shown in fig. ^ 
as a function of r^NLSp/''^ for rriNLSP = 100, 300 GeV and ITeV. Once both Onlsp^^ and 
'TiNLSP are determined from experiment, the upper bound on mLSP can be derived. Note 
that fig. P actually applies to both the axino and the gravitino LSP since it follows from 
eq. ( p.?!) . On the other hand, the bound is interesting only if f^NLSP^^ ^ 0.1, otherwise it 
reduces to a trivial condition mLSP < "znlsp- 

In fig. ^ we plot in the plane of varying rriNLSP and r^NLSP^^ where contours of Tr 
for a fixed ma = 0.01 GeV (left panel) and 1 GeV (right panel). In the left panel, where 
rria is tiny, for sufficiently large Tr thermal production dominates and the dependence on 
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Figure 4: Contours of the reheating temperature in the plane of tonlsp and JInlsp/i^ such that 
^lah'^ = flcDuh'^ = 0.104. The axino mass is assumed to be 0.01 GeV (left panel) and 1 GeV (right). 

f^NLSP^^ is very weak. Indeed, only at small rua and large J^nlsp^^ the NTP contribution 
oc (ma/mNLSp)^NLSP^^ starts playing some role. 

On the other hand, in the right panel of fig. Q the situation is somewhat more complex, 
with a characteristic turnover which depends on Tr. In order to explain it, in fig. |^ the 
TP and NTP contribution to the axino relic abundance are shown explicitly for one of the 
cases presented in fig. along with the corresponding Qnlsp/i^/IOO. As the neutralino 
mass = jtt-nlsp increases, the gluino mass also increases (via the gaugino unification 
mass relation). So long as niy^ < 130 GeV, the gluino decay contribution to Qj^h"^ dom- 
inates over the one from squark decays (with squark masses fixed in here at 1 TeV) and 
from scatterings (compare the left panel of fig. |^), but it decreases from left to right. 
This is because it is proportional to rrig times an integral over the Boltzmann suppression 
factor exp (— m^/r). On the other hand, the squark decay contribution, which scales as 
[m^ (1 — 0.05 log(mg/l TeV))]^, increases |l^. Overall, as increases, the TP part of 
JIq/i^ first decreases before increasing again. Thus the NTP part has to first increase and 
then decrease in order for the sum to remain constant at 0.104. 

In the upper left-hand corner of the right panel of fig. |^, for a given ratio of ?7T-a/mNLSP; 
^^NLSP^^ becomes too large for Qah'^ not to exceed 0.104, despite basically turning off the 
TP contribution (by reducing Tr). As Tr increases, the Boltzmann suppression factor 
plays a smaller role and the TP part increases, thus making the NTP part decrease. This 
is why, for the same value of ttinlsP) the corresponding value of I^nlsp^^ decreases with 
increasing Tr, in agreement with the right panel of fig. ^ Generally, for rria = 1 GeV the 
allowed values of Tr are much smaller than at smaller (left panel). 

In both figs. || and ^ we can see that a whole spectrum of values of J^nlsp^^ coming 
out from LHC determinations can be consistent with the axino relic abundance satisfying 
JIq/i^ ~ 0.1. They can range from much smaller to those much larger than the "WMAP 
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Choi, Roszkowski, Ruiz '07 

a LSP 
mg=lGeV - 
TR=110GeV ~ 




400 



Figure 5: The TP (dashed red) and NTP (dotted blue) contributions to the axino rehc abundance 
vs. mNLSP = for the case of ma = 1 GeV and Tr — 110 GeV from fig. ^. The CDM relic 
abundance has been set at its central value of 0.104. Also shown is fiNLSP^^/lOO (solid black) for 
this case. 



range" (IT), depending on the mass of the LSP. Moreover, f^NLSP^ does not have to 
take largest values in the NTP-dominated case, as one would initially expect. On the 
other hand, when mNLSP is much larger than ttt-lsp, NTP typically dominates, especially 
if i^NLSP^^ is also large. On the other hand, if /tt-nlsp is not much larger than mLSP and 
^NLSP^^ is also on the small side, then TP typically dominates. While this dependence is 
not rigorous, it may help one distinguishing between the TP and NTP-dominated regimes 
once mNLSP and J^nlsp^^ are experimentally determined. 

Finally we comment on the dependence on fa- Its value is limited to lie above some 
10^ GeV from cooling of red giants and below some 10^^ GeV from inflation if it took place 
prior to the Peccei-Quinn transition. Otherwise larger values are allowed |^0|. In our 
numerical examples we have set fa = 10^^ GeV. Since only TP of axinos depends on fa (as 
flj^h'^ oc 1//2) while NTP does not, in figs. | and | the regions dominated by TP will be 
affected accordingly. In other words, Tr oc f^/rtia, as stated above. 



4. Gravitino dark matter 

The case of gravitino LSP, with its relic abundance Q^h?' satisfying the range ( |1.1[ ), can 
be analyzed in an completely analogous way, except for two important aspects. 

The first has to do with Big Bang Nucleosynthesis (BBN). After the NLSPs freeze-out 
from thermal plasma and start decaying into the LSPs, energetic photons and hadronic 
jets are also produced. If this takes place during or after BBN, such particles may destroy 
successful predictions of standard BBN of light element abundances. In the axino LSP 
case, the NLSP decay occurs normally before BBN epoch begins |^, Thus the axino 
LSP is mostly free from BBN constraints. However, with gravitino couplings to ordinary 
matter being so much weaker, NLSP decays to gravitinos take place much later, varying 
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Figure 6: Left panel: Tr vs. ^NLSph'^ for Jtinlsp = 300 GeV and for mg = 0.01 GeV (solid blue) 
1 GeV (dashed red). The bands correspond to the upper and lower limits of dark matter 



and itIq 



density in (|l.l| ). Right panel: Tr, vs. for ^TssLSph — 100 (dashed blue), 0.1 (solid red) and 0.01 
(dotted black). To the right of the solid vertical line the gravitino is no longer the LSP. 



from 1 sec to 10^^ sec, depending on the masses of the two. Thus BBN usually provides 
severe constraints on the yield of the NLSP. 

The neutralino NLSP is almost excluded with > 1 GeV due to BBN j2|, |2|]. On 
the other hand, if the gravitino mass is less than ~ 1 GeV, the lifetime of the neutralino 
becomes smaller than about 1 sec and the window of neutralino NLSP opens up again |28|. 
In our numerical examples below we will not impose the constraint from BBN but will 
indicate in which case it is satisfied or not. 

The second important difference with the axino case was already discussed in sec. |2|. 
While is independent of the axino mass, in the gravitino case oc l/mg2. Thus 
rjj^/i^ oc rriaTji while Q^h^ oc T^/m^. In other words, if TP dominates, Q^h^ ~ 0.1, 
we find 



Tr cx mg. (4.1) 

This is reflected in the right panel of fig. ^ where we show Tr vs. for our nominal 
value of the NLSP mass of 300 GeV and for some fixed values of Onlsp^i^- So long as 
<^ ?7iNLSP find r^NLSP^^ is small enough, e.g. ^^nlsp^^ = 0.1 (denoted by a solid 
red curve), or less, the gravitino relic density is dominated by TP and Tr scales linearly 
with uIq so that Oj^/i^ ~ 0.1. On the other hand, when JInlsp^^ = 100 (dashed blue), 
NTP becomes rapidly dominant even with sub-GeV LSP mass. This can also be seen 
in the left panel of fig. |6| where we show Tr vs. I^nlsp^^ for the same NLSP mass and 
for mg = 0.01 GeV (solid blue) and mg = 1 GeV (dashed red). For the smaller value of 
uIq the TP contribution remains dominant up to much larger allowed values of J^nlsp^^ 
(horizontal part of the curves) before finally NTP takes over (vertical part). Note that, 
just like with the axino LSP case before, once we know ttj-nlsp aiid Onlsp^^j we can place 
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Figure 7: Maximum relieating temperature T^'^^ vs. NLSP relic density JInlsp^^ with gravitino 
DM for NLSP mass ttinlsp = 100 GeV (dashed red) and 300 GeV (sohd blue). 

an upper bound on the mass of the gravitino assumed to be the dominant component of 
CDM in the Universe, see fig. |3[ 

The turnover between the TP and NTP dominance allows one to derive a conservative 
upper bound T^^^ which, unlike for the axino CDM, even without knowing the gravitino 
mass. This can be seen from eqs. ( p.6D ~(p^) from which an expression for Tr as a quadratic 
function of mg can easily be found. Its maximum can then be expressed as 

for r^NLSP^^ > 0.05 

' (4-2) 

for ^NLSp/i < 0.05 

where we have used Qlsp^^ = 0.1. This is plotted in fig. |^ for two representative choices 
of TTT-NLSP = 100, 300 GeV, the former one being roughly the lowest value allowed by LEP. 

In fig. ^ we present, in analogy with fig. |^, contours of Tr for a fixed = 1 GeV 
(left panel) and 100 GeV (right panel). In both, at sufficiently large ^nlsp^'^ we find a 
similar turnover of the contours as in the right panel of fig. ^ for the axino case. Again, in 
order to elucidate the situation, in fig. ^ the TP and NTP contributions to the gravitino 
relic abundance are explicitly shown for one of the cases presented in fig. |^ along with 
the corresponding f^NLSP^^/lOO, in analogy with fig. ^ for the axino case. The (initially 
subdominant) TP part increases with increasing rUg (which grows with m^), thus the NTP 
part has to decrease in order for the sum to remain constant at 0.104. In the upper left- 
hand corner of both panels of fig. for a given ratio of m-g/w-NLSPi ^^NLSP^^ becomes 
too large and Q.^h? exceeds 0.104, despite basically turning off the TP contribution (by 
reducing Tr). 

In the gravitino LSP case it is easy to derive a simple expression for the value of the 
NLSP mass at which the turnover in I^nlsp^^ in fig- % takes place. By taking the NLSP to 
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Figure 8: Contours of the reheating temperature in the plane of ttinlsp and fiNLSp'i^ such that 
Vl^h? = ilcDM^^ = 0.104. The gravitino mass is assumed to be 1 GeV (left panel) and 100 GeV 
(right panel). 

Oo Choi, Roszkowski, Ruiz '07 
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Figure 9: The TP (dashed red) and NTP (dotted blue) contributions to the gravitino relic abun- 
dance vs. rriNLSP = '^^^ '^g ~ ^ ^"^"^ ~ ^^'^ from fig. ||. The CDM relic 
abundance has been set at its central value of 0.104. Also shown is riNLSp/i^/100 (solid black) for 
this case. 



be the bino (as we do in our numerical examples), plugging the expression for Q^h? from 



eq. ( p.6p into eq. ( |2.9D and making a simplifying assumption mg{^) ~ nig ~ Q.bm^^ one 
obtains an expression for Q^/i^ as a cubic function of m^. By putting its first derivative to 
zero one arrives at 



~ 54 



which adequately reproduces the positions of the peaks in fig. 
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The case of mg = 1 GeV is on the borderhne of being consistent with BBN, while the 
mass of 100 GeV is not (if the NLSP is the neutrahno). On the other hand, both cases 
can be allowed for some choices of MSSM parameters if the NLSP is the stau, see below. 
Note that in the gravitino case resulting values of Tr are typically significantly larger than 
for axinos. On the other hand, for sufficiently small < 1 GeV (roughly consistent 
with BBN constraints for NLSP neutralino) we find Tj^ ~ 10^~^ GeV, while for MeV-mass 
axinos Tr ~ 10^~^ GeV (compare the right panel of fig. which are of similar order. 

5. Stau NLSP 

Finally, we discuss the case of the lighter stau as the NLSP. Its couplings are of a sufficiently 
similar size to those of the neutralino that, for the same mass, its freeze-out temperature 
and yield are basically the same, hence is also its NTP contribution to either axino or 
gravitino relic abundance. For this reason, our numerical results presented in figs. for 
the neutralino NLSP case apply directly to the stau NLSP case.^ 

What is different is the impact of BBN constraints. Like with the neutralino, stau 
decays result in both electromagnetic and hadronic showers. Axino production from NLSP 
decays mostly takes place before the period of BBN, thus this case remains basically un- 
constrained. On the hand, the gravitino LSP case is again strongly constrained. Detailed 
analyses have shown that BBN constraints are somewhat weaker with stau NLSP than 
for the neutralinos |26, 27, 28|. As a result, the LSP gravitino mass in the range of tens, 
or even hundreds, of GeV, can remain allowed and consistent with Q.^h? in the "WMAP 
range" . Stau NLSP is also disallowed if its lifetime is longer than about 10^ sec, as otherwise 
it would form metastable bound states with nuclei which would affect BBN abundances 



of light elements, as recently pointed out in |31], see also [32, 33]. The stau lifetime is 
dominated by two-body decay to tau and gravitino which, neglecting tau mass, is given 
by [H 

rif^r + G) = 6.1 x 10^ sec ('-^\ ' (^^Y f 1 " ^ ■ (^-D 
^ ^ \ mf J V 100 GeV/ I m? ' ^ ^ 



In our numerical example in the right panel of fig. with mf = 300 GeV, where we 
have also taken = 477 GeV, the condition Tf > 10'^ sec implies < 2 GeV and 
^ 9 ^ 10^ GeV. Increasing nif to 1 TeV and to 1.5 TeV leads to < 40 GeV and 
Tr < 4 X 10^ GeV. 

The main point is that, with the neutralino NLSP, BBN constraints basically exclude 
gravitino CDM with the mass larger than about 1 GeV, while in the stau NLSP case much 
larger values can be allowed. 



6. Summary 

If the axino or the gravitino E-WIMP is the lightest superpartner and the dominant com- 
ponent of CDM in the Universe, then, under favorable conditions, it will be possible to 

*One caveat is that the gluino mass, which is important in TP, is not directly related anymore to the 
stau mass. 
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determine the reheating temperature after inflation from LHC measurements. Basicafly, 
one will need to measure the neutralino mass and other parameters determining its "relic 
abundance" - a non-trivial and highly exciting possibility in itself. The quantity may be 
found to be different from the "WMAP range" ( [1.1| ) , thus suggesting a non-standard CDM 
candidate. Additionally, one will need to know the mass of the gluino, either through a di- 
rect measurement, or via the gaugino uniflcation mass relation, since it plays an important 
role in E-WIMP production at high Tr. The neutralino is likely to be sufficiently long-lived 
to appear stable in LHC detectors, but in the early Universe it would have decayed into the 
axino or the gravitino LSP. Since each of them is extremely weakly interacting, it would 
be hopeless to look for them in usual WIMP dark matter searches. 

In the axino LSP case, if thermal production dominates, a consistency of its relic 
abundance with the "WMAP range" ( |1.1D implies Tr oc f'^/nia- If non-thermal production 
is dominant instead, then, for a known value of m-^^ the value of Q^^h? will imply an upper 
bound on ma/ fai which in turn will imply a lower bound on Tr. 

For the gravitino LSP, in the thermal production dominated case the analogous de- 
pendence is Tr oc mg. If non-thermal production dominates then the value of ^^^"^ 
imply an upper bound on m^^. For both axino and gravitino one can express this as 

max • r ^^CDM^l^ . 

"^LSP = mmjmNLSP, 7; ri^T-NLSpj, (6.1) 

which is independent of Tr. 

For both relics, one can also derive an upper bound on Tr, although in the axino 
case has to further assume that it is cold DM. On the other hand, for both axino and 
gravitino, it will be difficult to distinguish between the TP and NTP regimes but some 
information may be helpful. Typically, when is large and Onlsp^^ is also large then 
NTP is dominant, while if J^nlsp^^ is small then TP typically dominates. 

Similar conclusions as for the neutralino (for both E-WIMPs) apply if the NLSP is 
the stau, instead. Additionally, the very discovery of a charged, massive and (seemingly) 
stable state at the LHC will immediately imply that dark matter in the Universe is not 
made up of the usual neutralino. 

With so many similarities between both E-WIMPs, a natural question arises whether 
it would be possible to distinguish them in a collider experiment. Unfortunately, with 
neutralino NLSP this is unlikely to be possible. On the other hand, with the stau NLSP 
being electrically charged, enough of them could possibly be accumulated. By comparing 
their 3-body decays, it may be possible to only tell whether Nature has selected the axino 
or the gravitino as the lightest superpartner, but to actually also determine their mass with 
some reasonable precision |34, E^]. Needless to say, this would allow one to determine also 



the reheating temperature. 

In this work, we focused on demonstrating the principle of determining Tr at the 
LHC, and did not concern ourselves with uncertainties of experimental measurements. 
These are likely to be substantial at the LHC, the case we have mostly focused on here, 
unless one considers specific models A more detailed study will be required to assess 
their impact. An analogous study can be performed in the context of the planned Linear 
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Collider. We also implicitly assumed that only one relic dominates the CDM component 
of the Universe, which in principle does not have to be the case. For example, the axion 
and the axino could co-exist and contribute comparable fractions to the CDM relic density. 
(Note that this could be the case also in the "standard" case of the stable neutralino.) This 
would introduce an additional uncertainty and would lead to replacing the derived values 
of Tr, with an upper bound on the quantity. 
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